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Abstract 

Einstein's theory of Brownian motion is revisited in order to formulate generalized kinetic theory 
of anomalous diffusion. It is shown that if the assumptions of analyticity and the existence of the 
second moment of the displacement distribution are relaxed, the fractional derivative naturally 
appears in the diffusion equation. This is the first demonstration of the physical origin of the frac- 
tional derivative, in marked contrast to the usual phenomenological introduction of it. Furthermore, 
Einstein's approach is generalized to nonlinear kinetic theory to derive the porous-medium-type 
equation by the appropriate use of the escort distribution. 
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In his landmark 1905 paper Einstein has presented a new derivation of the diffusion 
equation. In this derivation, there appear two kinds of distributions. One is the probability, 
f(x,t)dx, of finding a Brownian particle in an interval [x,x + dx] at time t, and the other 
is the probability density, 0(A) for displacement, A, of the particle within a single discrete 
time step, r. The basic evolution equation is the following: 



/oo 
dAf(x + A, t)0(A) , 
-oo 



(1) 



where 0(A) satisfies the condition, 0(A) = 0(— A). It is clear that this equation is consistent 
with the normalization conditions on f(x,t) and 0(A). Assuming analyticity of f(x,t) in 
terms of both x and t and existence of the second moment of 0(A), Einstein derived the 
diffusion equation, df/dt = Dd 2 f/dx 2 , where the diffusion constant is given by 

A 2 



1 r°° A z 
D = - dA— 0(A) . 
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As well known, the mean displacement accordingly grows in time as 



(2) 
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where the over-bar denotes the expectation value with respect to f(x,t). It is important to 
note that Eq. (1) is different from today's familiar method of Green's function, f(x,t + r) = 
JZo dx'G{x, t + t\x', t)f{x f , t). 

Now, in many natural systems, we often observe diffusion processes, which do not follow 
the law in Eq. (3) but rather the spatial spread A oc t a with a different from 1/2. Such 
a phenomenon is referred to as anomalous diffusion, which is of extreme general interest 
in contemporary statistical mechanics Examples exhibiting anomalous diffusion 

are motion of tracer particles in turbulent flows ^j], charge transport in anomalous solids 

, porous glasses [10J, and 

subrecoil laser cooling [11]. To incorporate this exotic phenomenon with the framework of 
statistical mechanics, it seems necessary to generalize traditional kinetic theory. In fact, 
there are two different app roaches to this problem known in the literature as fractional 
H, Q and nonlinear 3, Q, Q Q, Q ones. To our knowledge, a 
connection between them has not been established yet on physical grounds. 

In this paper, we aim to derive and thereby unify the fractional and nonlinear theories 
of anomalous diffusion by appropriately generalizing Einstein's approach. A key point here 



is that the basic equation, Eq. (1), is an integral equation. This fact allows to relax the 
analyticity conditions on the relevant quantities. Another point is to include nonlinearity in 
Eq. (1) in conformity with the normalization condition, which may account for a nontrivial 
structure of the medium. 

Let us recall that Einstein assumed the existence of the second moment of the distribution, 
0(A), to derive the ordinary diffusion equation. This, however, does not hold for a large 
class of power-law distributions, like the Levy stable distributions, in particular. In this 
case, the expansion of f{x + A, t) in terms of A does not make sense. Instead, here we note 
that the right-hand side of Eq. (1) has the form of convolution of / and 0, which factorizes 
in the Fourier space: 

/(M + r)=/(M)#) , (4) 

where f(k,t) and (f)(k) are the characteristic functions of f(x,t) and 0(A), respectively, 
provided J 7 (g)(k) = g(k) = /f^, dy g(y)e tky . Now, suppose as an example of the distribution 
with the divergent second moment the Levy distribution, whose characteristic function is 
given by the stretched exponential form 

0(Jfe) = e - a W~' = l- a |fc|7 + ... ; (5) 

where a is a positive constant and the Levy index, 7, is in the range (0,2). 0(A) decays 
as a power law: 0(A) ~ |A| -1-7 . Substituting Eq. (5) into Eq. (4), performing the inverse 
Fourier transformation, and then expanding the left-hand side in terms of r, we obtain for 
the leading terms the following fractional diffusion equation: 

df <9 7 f 

-L = p* u J (Q) 

dx d\x\t ' w 

where the generalized diffusion constant, D*, is calculated to be 

D* = - . (7) 

T 

In deriving Eq. (6), we have used the the mathematical fact that the fractional derivative 
is defined through the Fourier transformation as 

^ (¥)<*) = -WTO • W 

The solution to this equation with the singular initial condition, f(x, 0) = S(x), is given by 

f(x,t) = / dke~ ikx - D * mi , (9) 

Z7T J -00 



which itself is the Levy distribution. From the scaling law, f(x,t) = t' 1 '^ f(x/t 1 ' 7 ), satisfied 
by the solution in Eq. (9), the spatial spread, A, (given by the half- width, for example) is 
seen to be, A oc t 1 / 7 , showing anomalous diffusion. It has to be emphasized that although 
fractional calculus has already been used for the discussions about Levy flight, the above 
derivation is the first to demonstrate its physical origin in a natural manner. 

Next, we proceed to show that Einstein's approach also provides a framework for diffusion 
phenomena in media with nontrivial structures with a unique viewpoint. In this case, we 
note that the structures can be taken care of by introducing nonlinearity in Eq. (1). Thus, 
we consider the following generalization: 

/oo 
dAF(f(x + A,t))<j)(A) , (10) 
-oo 



where F(f(x,t)) stands for the escort distribution |24 12 

F(f(xt))= $(/0M)) (11) 

[m)) JZodx'Hf(x>,t)) {LL) 

with $(y) a positive differentiate function. 

Let us discuss a special case when the medium has the following property. Suppose the 
total medium is divided into two uncorrelated subsystems, A and B, and their physical 
properties are identical to that of the total medium. In this case, the function, $, should 
satisfy the following condition: 

$(a&) = $(a)$(6) , (12) 

which implies that $ is a pure power-law function and correspondingly the escort distribution 
has the form 

^■'>> = r/w) ' (13) 

Now, assuming the existence of the second moment of 0(A), and then following Einstein's 
procedure, we arrive at a family of equations 

i*M = |j , (15) 

where c(t) = ^ oa dxf u {x 1 t) and D is given in Eq. (2). Except for the second term on 
the right-hand side as well as the time dependence in the generalized diffusion coefficient, 



£>t(t), E q. 



4) is the porous medium equation 19], which exhibits anomalous diffusion, 
A oc t~ (21 ] . The time dependence in the generalized diffusion coefficient may describe 
variations of the structure of the medium, seen from the diffusing particle. It should be 
calculated by using the solution of Eq. (14) in a self-consistent manner. The second term 
on the right hand side can be large due to the factor, 1/r. However, the spatial integral 
of this term identically vanishes because of the normalization conditions. Therefore, it is 
regarded as a fluctuating contribution. Unfortunately, this novel equation does not seem to 
be analytically tractable, and heavy numerical analysis is yet to be carried out to examine 
the detailed diffusion properties. 

In conclusion, we have revisited Einstein's 1905 theory of Brownian motion in a modern 
context of anomalous diffusion. We have discussed both fractional and nonlinear general- 
izations of ordinary kinetic theory. We have shown how naturally fractional calculus enters 
into the generalized diffusion equation if the assumptions of analyticity and the existence 
of the second moment are relaxed. To our knowledge, this is the first demonstration of the 
physical origin of the fractional derivative, in marked contrast to the usual phenomenolog- 
ical introduction of it. We have also generalized Einstein's approach to nonlinear kinetic 
theory. We have seen that the porous-medium-type equation can be derived if the escort 
distribution is appropriately employed. 

It is a pleasure to see how Einstein's century old idea can still provide a new insight into 
a modern aspect of Brownian motion. 
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